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Abstract 

In generally covariant metric gravity theories with tensor matter 
fields, the initial value constraint equations, unlike in general rela- 
tivity, are in general not just the Ofj, components of the metric field 
equation. This happens because higher derivatives can occur in the 
matter stress tensor. A universal form for these constraints is derived 
here from a generalized Bianchi identity that includes matter fields. As 
an application, the constraints for Einstein-aether theory are found. 

1 Introduction 

A field theory with diffeomorphism symmetry can not have deterministic 
field equations. Although Einstein briefly thought this meant that a phys- 
ical theory could not be diffeomorphism invariant, he soon realized that 
the observables may be deterministically evolved, even though the fields are 
notpQ). The solutions must involve four free functions (in four spacetime 
dimensions), which correspond to the freedom to apply an arbitrary time 
dependent diffeomorphism to the fields. This means that out of the col- 
lection of Euler-Lagrange equations obtained from varying the fields in the 
action, there must be four combinations that are not independent of the 
others. These are initial value constraint equations, and the remaining field 
equations imply that if the constraints hold initially, they are automatically 
preserved in time. 

This is all very familiar to anyone who has worked with general rela- 
tivity, but there is a twist if the matter energy-momentum tensor involves 
second time derivatives of the "matter" fields, as happens if the matter ac- 
tion involves covariant derivatives of tensor fields. In that case, the initial 
value constraints are not just the 0/i components of the metric field equa- 
tions; instead, one must add to those a certain combination of the matter 
field equations. An example is Einstein-aether theory, in which the metric 
is coupled to a unit timelike vector field. 
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In this article the form of the constraint equations in the presence of 
tensor matter is derived from a generalization of the Bianchi identity that 
includes the matter fields. It is then applied to three examples that illustrate 
different features. In particular, the constraints for Einstein-aether theory 
are obtained. After writing this up I learned that the result is (of course) 
not new. A paper by Bergmann[2] from 1949 obtains the same result, and 
it probably goes back even further. It was recently derived by Seifert and 
Wald[3], and applied to Einstein-aether theory in [3]. However, it appears 
not to be common knowledge today. Since theories with tensor matter are 
currently of interest, it seems worth publishing this brief note, which aims to 
explain the result and make it easily accessible using low brow, coordinate 
based methods. 



2 General relativity 

To begin with let us review how it works in general relativity. 

2.1 Vacuum 

The vacuum Einstein equation, 

G a/3 = 0, (1) 

involves 10 independent components, matching the 10 independent compo- 
nents of the spacetime metric. However 4 of the 10 equations packaged into 
([1]) must not actually be evolution equations. In fact, for any coordinate 
system x a = (x°,x l ) (i = 1,2,3), the 0/3 components of (pQ) involve no sec- 
ond derivatives with respect to x°. (I will refer to x° as the "time" in the 
following, although in fact it can be any coordinate.) A simple way to see 
thiai] is from the Bianchi identity, 

V a G a/3 = 0, (2) 

which in component form reads 

d G 013 + diG ifi + r a aa G a P + TP afT G aa = 0. (3) 

lr The argument is given this way in Weinberg's 1972 book[5]- It also appears in the 
1975 english edition of Landau and Lifshitz's vol. 2[6], which is translated from the 1973 
russian 6th edition, but not in the 1971 english edition. For a chronicle of the early history 
see Stachel0. 
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were to have second time derivatives, then the first term in ([3]) would 
have third time derivatives, unlike any other term. This would be inconsis- 
tent with ([3]) being true as an identity. Therefore we can infer that G°P has 
no second time derivatives. The equations 

G 0/s = (4) 

are therefore not evolution equations but rather constraints on initial data. 
Moreover, the Bianchi identity is a first order differential equation doG @ = 
. . . , which implies that if these constraints hold at some initial time, and 
the rest of the field equations G lJ = hold at all times, then the constraints 
hold for all times. 

It is no accident that the Bianchi identity is so intimately tied to the 
existence and preservation of the constraint equations. That identity can 
be viewed as a direct consequence of the diffeomorphism invariance of the 
Einstein-Hilbert action, S = f R. To reduce clutter, the metric volume 
element y/—g d 4 x is implicit here and in the following. (When the metric 
is varied, yj—g must certainly be varied as well.) Under a metric variation 
8g a p the action varies a^l 

SS = - J G afi 5g aP (5) 

If the metric variation arises from an infinitesimal diffeomorphism generated 
by a vector field we have 

%a/3 = = 2V( a ^) (6) 

where C denotes the Lie derivative, and the round brackets on indices denote 
symmetrization. Substituting (JSJ) in (JSJ), integrating by parts, and assuming 
vanishes at the boundaries, we find that a diffeomorphism generates the 
variation 

S^S = J (2V a G°%. (7) 

Since the Einstein-Hilbert action is a spacetime scalar, 5^S must vanish for 
all which implies the Bianchi identity (J2J). That is, the action is au- 
tomatically stationary under those variations of the metric that arise from 
diffeomorphisms, so the corresponding Euler-Lagrange equations do not im- 
pose any conditions. The identity implied by this symmetry indicates which 
four of the field equations are constraints rather than evolution equations, 
and ensures that those constraints are preserved in time. 



2 We use the standard notation: indices are lowered and raised with the metric and 
inverse metric, R = g a ^ R a p is the Ricci scalar, R a p — R" aa p is the Ricci tensor, R" a p-y 
is the Riemann tensor, and G a p = R a p — ^Rg a p is the Einstein tensor. 
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2.2 Matter 



Now let us consider what happens when matter is coupled to the metric. 
The field equations then take the form 

G afS = 8vrT a/3 (8) 

(in units with c = G = 1), where T a ^ is the matter stress tensor. If the stress 
tensor involves no higher than first derivatives and the matter field equations 
are of second order in derivatives, then the equations G 0/3 = 8ttT ^ are again 
initial value constraint equations. However, if the matter action involves 
the covariant derivative, such as for example in a term like (V a n /3 )(V Q n /3 ) 
(where u 13 is a vector field), then variation of the metric will produce (among 
other things) terms involving the second derivatives of the metric and matter 
fields, so T a P will generally involve second time derivatives. Then G ^ = 
8ttT ^ will not be constraint equations. 

The difference from the vacuum case might seem puzzling since, when the 
matter satisfies its equation of motion, its stress tensor is divergence-free, 
V a T a/3 = 0. Together with the Bianchi identity this yields the equation 
V a (G a P — 87rT a P) = 0, and it is tempting to argue by analogy with the 
Bianchi identity case that validity of this equation implies that G ^ — 8ttT ^ 
can have no second time derivatives. This would be incorrect however, since 
the equation uses the equations of motion of the matter field, so does not 
hold identically for all matter fields. In order to discover the true constraint 
equations in theories like this, what we need to do is find the analog of the 
Bianchi identity, that is, the identity that follows from the diffeomorphism 
invar iance of the action. 

3 Any metric theory with tensor matter 

Consider then a diffeomorphism invariant action S[g, ip] that is a functional 
of the metric g a p and some tensor matter field of any rank, ij)a.'.'.- For 
notational simplicity I will assume ip is a (1, 1) tensor in what follows but 
it should be obvious how to modify the equations to allow for the general 
case. A diffeomorphism induces the variation ([6]) for the metric and 

5^ a = = c v 7 v$ + < v a c - C v 7 ^ (9) 

for the tensor. Actually the Lie derivative is independent of the metric, 
so all the Christoffel symbols in ([9]) cancel out, but to maintain manifest 
covariance and to avoid the need to explicitly write the yj— g factor it is 
more convenient to use covariant derivatives. 
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The variation of the action under a diffeomorphism is 

<%S = J E a ? S i9afi + H$ S^ a , (10) 

where 

= H t = (11) 

E a @ = is the metric field equation, and Hp = is the matter field equation. 
Now we insert the diffeo variations ([6]) and ([9]) into (I10p . integrate by parts 
on all derivatives of £ Q , and redefine some dummy indices to obtain 

J [V a (-2E% - Hy} + Hj^) - H«V^l] f. (12) 

Since the action is diffeomorphism invariant, this variation must vanish for 
all vector fields so we infer an identity, 

V a {-ZE$ - H^J + HjW = H« V^l (13) 

This generalizes the Bianchi identity ([2]). Note that, unlike the latter, it is 
not just the statement that a certain tensor is divergence-free, but rather 
that the divergence is equal to some other quantity. It includes the matter 
fields and holds off-shell, i.e. for all field configurations. 

Now to identify the initial value constraint equations, we expand the co- 
variant divergence on the left hand side of (|13p in terms of partial derivatives 
and Christoffel symbols as in (J3]). The term d (-2£° _ #0^,7 + m^°) will 
have one higher time derivative order than any other term in (|13p . Since 
(I13p holds identically for all fields, we infer that the quantity 

+ (14) 

has time derivative order that is one less than the highest order in the field 
equations E a @ = and H% = 0. Also, it is a combination of components 
of the field equations, so it vanishes when the field equations hold. We thus 
arrive at the conclusion that the initial value constraints are given by 

Cp = -2El - HfyJ + = 0. (15) 

For a general rank tensor matter field ip[" there would be one term in Cp 
for each index. Each covariant matter index contributes —H'.'.'P'"ij}\[[p___, and 
each contravariant index contributes iT'.'^...^.'.'.' "'. For multiple matter fields, 
there would be similar terms for each matter field. The generalized Bianchi 
identity (|13|) implies that the constraints are automatically preserved in time 
if they hold at one time and the matter field equations and other components 
of the metric field equations hold at all times. 
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4 Examples 



We now apply this result to three examples. The first example is Einstein- 
Maxwell theory. The matter field is a covariant vector potential A a , and the 
matter Lagrangian is — 4 F a ^F a p, with F af} = daAp - dpA a . The Maxwell 
stress tensor has no second derivative terms, so the quantities E 0/3 = —G°^ + 
8ttT°P involve no second time derivatives. Nevertheless, it is also true that 
Cg in (|15p has no second time derivatives. The extra term in Cg in this 
case is —H°Ag, so it must be that also H° has no second time derivatives. 
Indeed, H° = V a F a0 , and the antisymmetry of F a @ implies that the a 
index must be spatial. This field equation corresponds to V • E = 0, the 
initial value constraint for Maxwell's equations. 

For the case of Einstein-aether theory[8j the tensor matter is the aether 
vector u a , and the Lagrangian contains Christoffel symbols, for instance in 
the term (V a ug)(y a u^), so the stress tensor has second time derivatives. 
The constraint equations are 

C |ther _ _ 2E + u Hp = Q) (16) 

where Hg = is the aether field equation, which contains second time deriva- 
tives that cancel those in E® when the two are added with the coefficients 
prescribed in (I16p . A special feature of this theory is that the field equation 
has no component in the direction of vP , i.e. HgvP = 0. That is because 
u@ is constrained to be a unit timelike vector, hence all of its variations 
are orthogonal to it, ugbvP = oH Thus the 0-u component of the Einstein 
equation, E^u^ = 0, is an initial value constraint without any additional 
terms involving the aether field equation. In particular, on a hypersurface 
to which the aether is orthogonal, the uu component of the aether stress 
tensor has no second derivatives in the u direction. 

The final example is Hofava gravity [9j 110] . I will first describe it in the 
original 3 + 1 formulation with a preferred time, and then explain how the 
constraints work in its fully covariant formulation. In the original form, this 
theory has a fixed time function, defined up to reparameterizations, that 
determines a foliation of spacetime. The remaining symmetry is the time 
dependent spatial diffeomorphisms on the constant time surfaces. The field 
variables are the spatial metric, the shift vector, and the lapse scalar. I will 
discuss only the version of this theory in which the lapse iV is allowed to de- 
pend on the spatial coordinates, and in which the action contains all terms 
up to two derivatives, including[ll] a term proportional to (diN)(d 1 N)/N 2 

3 If the unit constraint is imposed by variation of a Lagrange multiplier, this statement 
holds only after the multiplier is solved for and inserted back into Hp. 
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(which is invariant under time dependent rescaling of N since <9j denotes 
the spatial gradient). This theory has three first class constraints that gen- 
erate spatial diffeomorphisms in the Hamiltonian formulation, and a fourth 
constraint that is second class, because the theory lacks surface deformation 
symmetry except for global time reparameterization [12 \ [T^ j 1 14] Fl This fourth 
constraint is the field equation for the lapse, and is an elliptic equation in 
spatial derivatives, so determines the lapse in terms of the other variables. 

In the covariant formulation |1 1 \ [T5] . the field variables are the spacetime 
metric and a scalar field T, with a symmetry under reparameterizations of 
T — > f(T). The two derivative action is just that of Einstein-aether theory, 
with aether vector u a = Nd a T and N = (g a/3 d a T dpT)^ 1 / 2 . Because the 
"matter" field T is a scalar, the generalized Bianchi identity (|13j) takes the 
form 



where H = -V Q [N(Sf - u a u )5S/5u p ] is the T field equation^. The 
covariant derivative of u a appears in the action so, although the matter is 
a scalar field, the aether part of Ea has derivatives of order two on g a p and 
three on T, while H has derivatives of one higher order on each field variable. 
It is not until the gauge x° = T is chosen that the time derivative order of the 
theory drops to two: the field variable T disappears from the equations, and 
the derivative orders of E^ and H remain two and three respectively, but the 
third derivative in the H field equation is always in a spatial direction. Since 
both field equations are now second order in time derivatives, the gauge- 
fixed identity (|17|) implies that all four quantities Eg have no second time 
derivatives. (The right hand side of (|1T|) has no component orthogonal to u a 
anyway, so even if H had third time derivatives, one could have still inferred 
that the time-space components Ef have no second time derivatives.) That 
is, the constraints are of the usual form, with no need to add a term involving 
the matter field equation. 

5 Discussion 

The constraint equations are essential input for calculation of evolution from 
initial data. That evolution could be in a timelike direction, as in ordi- 
nary time evolution, or a spacelike direction, as e.g. in the case of finding 
static or stationary solutions to field equations. In particular, the constraint 
equations for Einstein-aether theory on a timelike surface of constant ra- 
dial coordinate were used in [16] to find static, spherically symmetric black 

4 On a compact spatial manifold there is one global first class surface deformation 
constraint that generates time reparameterizations [13] ■ 




(17) 
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hole solutions to that theory. The constraints on a spacelike hypersurface in 
spherical symmetry were used in [4] to study perturbative stability of solu- 
tions in various modified gravity theories including Einstein-aether theory. 
The constraints for a spherically symmetric, spacelike hypersurface to which 
the aether is orthogonal were worked out for Einstein-aether theory in |17j . 
and used in a numerical simulation of collapse of matter to form a black hole. 
An interesting application of the generalization beyond spherical symmetry 
would be to numerical simulation of axisymmetric black hole collapse, which 
could yield at late times a stationary rotating black hole solution. 
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